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Solutions to the Atmospheric Neutrino Problem 
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In this talk I review the present status of the atmospheric neutrino anomaly and discuss some solutions that 
have been presented in the literature to solve this problem. In particular I review the "standard" solution in 
terms of neutrino oscillations as well as alternative scenarios such as the possibility of flavour changing neutrino 
interactions with the Earth and neutrino decay. 



1. Introduction 

Neutrinos produced as decay products in 
hadronic showers from cosmic ray collisions with 
nuclei in the upper atmosphere have been ob- 
served by several detectors ^^- Although 
the absolute fluxes of atmospheric neutrinos are 
largely uncertain, the expected ratio (n/e) of the 
muon neutrino flux {v^ + J'/i) over the electron 
neutrino flux (i>e + i^e) is robust, since it largely 
cancels out the uncertainties associated with the 
absolute flux. In fact, this ratio has been calcu- 
lated gl with an uncertainty of less than 5% over 
energies varying from 0.1 GeV to 100 GeV. On 
this resides our confidence in the long-standing 
atmospheric neutrino anomaly. 

Super-Kamiokande high statistics observa- 
tions H] indicate that the deficit in the total ratio 
R(/x/e) is due to the number of neutrinos arriv- 
ing in the detector at large zenith angles. The 
e-like events do not present any compelling ev- 
idence of a zenith-angle dependent suppression 
while the ^-like event rates are substantially sup- 
pressed at large zenith angles. The v^ —^ Vr as 
well as the z^^ -^ Vg Pp0| oscillation hypoth- 
esis provides a very good explanation for this 
smaller-than-expected ratio, which is also sim- 
ple and well-motivated theoretically. This led 
the Super-Kamiokande Collaboration to conclude 
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that their data provide good evidence for neu- 
trino oscillations and neutrino masses [^. How- 
ever, alternative explanations to the atmospheric 
neutrino data have been proposed in the litera- 
ture including the possibility of flavour changing 
(PC) neutrino interactions in matter ||l2| and neu- 
trino deca y [13| , the violation of relativity prin- 
ciples |lj,f5| or the violation of GPT symmetry 

El]- 

In this talk I will review the "standard" so- 
lution in terms of neutrino oscillations as well 
as alternative scenarios such as the possibility of 
flavour changing neutrino interactions with the 
Earth and neutrino decay. 

2. Atmospheric Neutrino Induced Events 
at Underground Experiments 

Atmospheric neutrinos can be detected in un- 
derground detectors by direct observation of their 
charged current interaction inside the detector. 
These are the so called contained events. Con- 
tained events can be further classified into fully 
contained events when the produced charged lep- 
ton (either electron or muon) in the neutrino in- 
teraction does not escape the detector, and par- 
tially contained muons when the produced muon 
exits the detector. Super-Kamiokande has di- 
vided their contained data sample into sub-GeV 
events with visible energy below 1.2 GeV and 
multi-GeV above such cutoff. On average, sub- 
GeV events arise from neutrinos of several hun- 
dreds of MeV while multi-GeV events are origi- 



nated by neutrinos with energies of the order of 
several GeV. Higher energy muon neutrinos and 
antineutrinos can also be detected indirectly by 
observing the muons produced in their charged 
current interactions in the vicinity of the detector. 
These are the so called upgoing muons. Should 
the muon stop inside the detector, it will be clas- 
sified as a "stopping" muon, (which arises from 
neutrinos of energies around ten GeV) while if 
the muon track crosses the full detector the event 
is classified as a "through-going" muon which is 
originated by neutrinos with energies of the order 
of hundred GeV.. 

Given certain neutrino conversion mechanism, 
the expected number of /x-like and e-like con- 
tained events, Na, a = fi,e can be computed as: 



where 



Ne^ 



Ne = Nee + N,, 



(1) 
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a/3 



TltT 



d2$„ 



dE^d{cos9^) 



Ka{h, cos9^, E^) 



da 



Po.f3-pp^e{Efj)dE^dEi3d{cos9^)dh (2) 
dEp 

and Pa/3 is the conversion probability of i/q —f Vfj 
for given values of Ei,,cos9„ and h, i.e.. Pap = 
P{va -^ Pp; E,y, cos 9v,h). In the Standard Model 
(SM) , the only non-zero elements are the diagonal 
ones, i.e. Paa = 1 for all a. Here Ut is the number 
of targets, T is the experiment's running time, E^, 
is the neutrino energy and $q. is the flux of at- 
mospheric neutrinos of type a = ii,e; Ep is the 
final charged lepton energy and s{Ep) is the de- 
tection efficiency for such charged lepton; a is the 
neutrino-nucleon interaction cross section, and 9^ 
is the angle between the vertical direction and 
the incoming neutrinos (cos 9^=1 corresponds to 
the down-coming neutrinos). In Eq. (y), /i is the 
slant distance from the production point to the 
sea level for a-type neutrinos with energy E^, and 
zenith angle 9^. Finally, Ka is the slant distance 
distribution which is normalized to one |17]| . 

The neutrino fluxes, in particular in the sub- 
GeV range, depend on the solar activity. In order 
to take this fact into account we use in Eq. (||) a 
linear combination of atmospheric neutrino fluxes 
^""^^ and (j)™*" which correspond to the most ac- 
tive Sun (solar maximum) and quiet Sun (solar 



minimum), respectively, with different weights, 
depending on the running period of each experi- 
ment ]9|. 

Upgoing muon data are usually presented in 
the form of measured muon fluxes. We ob- 
tain the effective muon fluxes for both stop- 
ping and through-going muons by convoluting the 
probabilities with the corresponding muon fluxes 
produced by the neutrino interactions with the 
Earth. We include the muon energy loss during 
propagation both in the rock and in the detec- 
tor according to [l^Jl^ and we take into account 
also the effective detector area for both types of 
events, stopping and through-going. Schemati- 
cally 
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where R(Ef^o, Ef^) is the muon range function 
which accounts for the muon energy loss during 



propagation. A{L, 



AsiE^,9) + AriE^ 



is the projected detector area for internal path- 
lengths longer than L. As and At are the cor- 
responding areas for stopping and through-going 
muon trajectories. For Super-Kamiokande we 
compute these effective areas using the simple ge- 
ometrical picture given in Ref. |20|] . 

Following Ref. |y| we explicitly verify in our 
present reanalysis the agreement of our predic- 
tions with the experimental Monte Carlo predic- 
tions, leading to a good confidence in the reliabil- 
ity of our results. 

3. Conversion Probabilities 

For definiteness I assume a two-flavour sce- 
nario. For the oscillation case one must solve the 
Schroedinger evolution equation of the v^, — vx 
(where X — e,T or s sterile) system in the mat- 
ter background for neutrinos 
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Here Gp is the Fermi constant, p is the mat- 
ter density in the Earth, M is the nucleon mass, 
and Ye (Yn) is the electron (neutron) fraction. I 



define Am = 



in such a way that if 



Am'^ > (Am^ < 0) the neutrino with largest 
muon-Hke component is heavier (hghter) than the 
one with largest X-like component. For anti- 
neutrinos the signs of potentials Vx should be 
reversed. In our calculations we have used the 
approximate analytic expression for the matter 
density profile in the Earth obtained in Ref. [gl[ . 
In order to obtain the oscillation probabilities Pa/s 
we have made a numerical integration of the evo- 
lution equation. The probabilities for neutrinos 
and anti-neutrinos are different because the re- 
versal of sign of matter potential. Notice that for 



the i^n 



case there is no matter effect and 



the probability takes the well known form 



P^^ ^ 1 -sin^26'sin2( 



Am?L 
2E 



)• 



(12) 



For the u^^ —^ Vg case there are two possibilities 
depending on the sign of Arn^. For Aw? > the 
matter effects enhance neutrino oscillations while 
depress anti-neutrino oscillations, whereas for the 
other sign {Am^ < 0) the opposite holds. In what 
follows 1 will not consider the possibility of oscil- 
lation into electron neutrinos as it is known to be 
ruled out by the negative results of the reactor 
experiment Chooz p^ . 

I am also going to consider the possibility of 
FC-neutrino interactions of massless neutrinos 
which can also induce ly,, -^ v-r transitions [tl2|. 



In our phenomenological approach we have as- 
sumed that the evolution equations which de- 
scribe i/^ — > Vt transitions in matter may be writ- 
ten as Eq.(H) with 

i?,, - Hr = V2 Gpelnfir) (13) 

H^r = V2 Gpe.nfir) (14) 

where \/2GFnf{r)€i, is the i^f^ + f ^ i^t + f for- 
ward scattering amplitude and y/2GFnf{r)e'^ is 
the difference between the Vr — f and v^ — J elastic 
forward scattering amplitudes, with nf{r) being 
the number density of the fcrmions which induce 
such processes. The parameters e^/ and e^ contain 
the information about FC-neutrino interactions, 
for details I refer to Ref. [O. In order to obtain 
the oscillation probabilities Pap we have made a 
numerical integration of the evolution equation. 
For the sake of illustration I show here the so- 
lution in the approximation of constant matter 
density. The conversion probability in this case 
is 



Puu = l- 



4e, 



4e2 



.'2 



HIvL), 



(15) 



where 77 = JAel + e^^V2G_Fn/. 

Finally when discussing the possibility of neu- 
trino decay in the Vf^ Vr system one must include 
also the effect of the neutrino instability when 
solving the time evolution of the mass eigenstates. 
In this case the evolution equation can be solved 
analytically and the muon neutrino survival prob- 
ability is given by: 

TO2 L^ 

~E~ 



Pi_t^, = sin' 



cos y exp 



Pa 



(16) 



with 



1 2 ( Arn^L 

= - sm 26* cos I ^^ ) exp 



m2 L\ 
~^27 J 



where r is the neutrino lifetime. For clarity, I 
isolated the term PAm^ i explicitly dependent on 
the mass difference, Arri^ = ra\ — rri^. In the 
following I will neglect this term which averages 
out in the limit mi ^ TO2 ^ E. The factor ni2/E 
in the exponential is just the 7 Lorentz factor. 

By comparing Eqs.(|lJ), (p^, and ( [iq ) one can 
see that the different mechanisms lead to differ- 
ent energy dependence of the survival probability. 



While oscillations give a muon survival probabil- 
ity slowly growing with the energy, neutrino de- 
cay leads to a much faster growth. On the other 
hand, flavour changing neutrino interactions lead 
to an energy independent survival probability. As 
pointed out in the previous section, data on con- 
tained events and upgoing muons scan a range 
of neutrino energies over three orders of magni- 
tude from hundreds of MeV to hundreds of GeV. 
This is of great importance when discriminating 
between the different scenarios as pointed out in 
Refs.fE3lJ2l. 



4. Atmospheric Neutrino Data Fits 

Here I describe our fitting method to determine 
the atmospheric neutrino conversion parameters 
for the various possible conversion channels. In 
doing so we have relied on the separate use of 
the event numbers paying attention to the cor- 
relations between the errors in the muon predic- 
tions and electron predictions as well as the cor- 
relations among the errors of the different energy 
data samples. 

The steps required in order to generate the al- 
lowed regions of oscillation parameters were given 
in Ref. ||. Following Ref. |,|| we defined the 
X^ as 






(iVf-iVf).(ai+4J7i.(iVf-Aff),(I7) 



where / and J stand for any combination of the 
experimental data set and event-type considered, 
i.e, / = [A,a) and J = (i?,/3) where, A,B 
stands for the different experiments or different 
data samples in a given experiment. a,/3 = e,^. 
In Eq. ( p7| ) Nf^ is the predicted number of events 
calculated from Eq. (^ whereas Nf'^ is the num- 
ber of observed events. In Eq. ( p^ ) cr^^ and cr^,^ 
are the error matrices containing the experimen- 
tal and theoretical errors respectively. They can 
be written as 



'ij 



(T^{A)p^p{A,B)(jp{B), 



(18) 



where pap{A,B) stands for the correlation be- 
tween the a-like events in the A-type experiment 
and /3-like events in B-type experiment, whereas 
aa{A) and cp^B) are the errors for the number 



of a and /3-like events in A and B experiments, 
respectively. The dimension of the error matrix 
varies depending on the combination of experi- 
ments included in the analysis. 

We have computed paf3{A,B) as in Ref. |2^ . 
A detailed discussion of the errors and correla- 
tions used in our analysis for the contained events 
can be found in Ref. [g[. In our present analysis, 
we have also included the data on stopping and 
through-going muons. We have conservatively as- 
cribed a 20% uncertainty to the absolute neu- 
trino flux, in order to generously account for the 
spread of predictions in different neutrino flux cal- 
culations. Other important source of theoretical 
uncertainty arises from the neutrino interaction 
cross section which at Super-Kamiokande ranges 
from 10-15 %. We allow a 5% variation in the ra- 
tio between muon events in different energy sam- 
ples. We further introduce a 10% theoretical error 
in the ratio of electron-type to muon-type events 
of the different samples. Uncertainties in the ratio 
between different angular bins are treated, simi- 
larly to Ref. Ig^ . With our definitions we obtain, 
for instance, Xsm =122/(35 d.o.f) which means 
that the SM has a CL of 10"" ! 

Next we minimize the x^ function in Eq. (|l7| ) 
and determine the allowed region in the parame- 
ter space for certain conversion mechanism and 
for a given confidence level, defined as, x^ — 
xLin + 4-61 (9.21) for 90 (99) % CL. 

In Table m I show the minimum value of x^ and 
the best fit point for several conversion mecha- 
nisms and for the different data sets. 

The results of our x^ fit of the Super- 
Kamiokande contained and upgoing atmospheric 
neutrino data in the framework of neutrino oscil- 
lations are given in Fig. (m). In this figure I give 
the allowed region of oscillation parameters at 90 
and 99 % CL. 

One can notice that matter effects lead to dif- 
ferences between the allowed regions in the dif- 
ferent channels. In the case of v^ — *■ I's with 
Arn^ > matter effects enhance the oscillations 
for neutrinos and therefore smaller values of the 
vacuum mixing angle would lead to larger con- 
version probabilities and the regions are therefore 
larger as compared to the vacuum case v^ ^ Vt- 
In the case of Vu — > Vs with A-m? < the enhance- 
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Figure 1. Allowed regions of the oscillation parameters for the different Super-Kamiokande data samples 
and oscillation channels as labelled in the figure. 



ment occurs only for anti-neutrinos while in this 
case the effect of matter suppresses the conversion 
in i^^'s. Since the yield of atmospheric neutri- 
nos is bigger than that of anti-neutrinos, clearly 
the matter effect suppresses the overall conver- 
sion probability. Therefore one needs in this case 
a larger value of the vacuum mixing angle, as can 
be seen by comparing the regions in the second 
row with the corresponding ones in the first and 
third row in Fig. (|l|). 

Notice that in all channels where matter ef- 
fects play a role the range of acceptable Ani^ is 
slightly shifted towards larger values, when com- 



pared with the Vu 



case. This follows from 



looking at the relation between mixing in vacuo 
and in matter. In fact, away from the resonance 
region, independently of the sign of the matter 
potential, there is a suppression of the mixing in- 
side the Earth. As a result, there is a lower cut 
in the allowed Am^ value, and it lies higher than 
what is obtained in the data fit for the v^ -^ v^ 
channel. 

Concerning the quality of the fits we see in ta- 
ble 1^ that the best fit to the full sample is ob- 
tained for the v^^ —> Vj- channel although from 
the global analysis oscillations into sterile neutri- 
nos cannot be ruled out. One can also observe 
an improvement in the quality of the fits to the 
contained events as compared to previous analy- 
sis performed with lower statistics [g[ . These fea- 
tures can be easily understood by looking at the 
predicted zenith angle distribution of the differ- 
ent event type for the various oscillation channels 
which 1 show in Fig. (0) for contained events and 
Fig. (ph for upgoing muons. In Fig. (0) we ob- 
serve a perfect agreement between the observed 
distributions of e-like events and the predictions 
in the SM. This has lead to an improvement of 
the quality of the description for any conversion 
mechanism that only involves muons. Also in 
Fig. (0) we can observe that due to matter ef- 
fects the distribution for upgoing muons in the 
case of Vf^ — > Vs are flatter than for v^ -^ Vt 
ll^,^. Data show a somehow steeper angu- 
lar dependence which can be better described by 
v^ —f Vr- This leads to the better quality of the 
global fit in this channel. Pushing further this 
feature Super-Kamiokande collaboration has pre- 



sented a preliminary partial analysis of the angu- 
lar dependence of the through-going muon data 
in combination with the up-down asymmetry of 
partially contained events which seems to exclude 



the possibility v^ 



at the 2-(T level 



5. Alternative Scenarios 

As we have seen in the previous section the 
oscillation hypothesis provides a very good ex- 
planation to the atmospheric neutrino data, and 
it is also simple and well-motivated theoretically. 
However, alternative explanations to the atmo- 
spheric neutrino data have been proposed in the 
literature. In this section I concentrate on the 
present status of two possible "exotic" scenarios: 
FC-neutrino interactions in the Earth matter |12] 
and neutrino decay p3| . 

In Fig. (H) I show the allowed regions in the 
parameter space for the decay mechanism v^ —^ 
VtX for Super-Kamiokande sub-GeV and multi- 
GeV events separately as well as the combined 
contained events. Partial fits to the two events 
samples lead to rather good description as can 
be seen in Table [2]. However the description of 



Table 2 

Minimum value of y^ and the best fit point to 

contained events for the decay Vn —^ VrX . 



Experiment 
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Amzn 

m/r (10- 
sin^e 


-5 GcV/Km) 


1.4 
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0.05 


Super-Kam 
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y' ■ 
m/r (10- 


-^ GcV/Km) 


5.9 
31. 
0.1 


Super-Kam 
Contained 


Arnin 

m/r (10- 

sin^ e 


'^ GeV/Km) 


21.9 

9.3 

0.07 



the global contained event sample is considerably 
worse than in the case of oscillations. This arises 
from the stronger energy dependence of the sur- 
vival probability while the contained data both 
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Figure 2. Angular distribution for Super- 
Kamiokande electron-like and muon- like sub- 
GeV and multi-GeV events together with our pre- 
diction in the absence of oscillation as well as the 
prediction for the best fit point to the contained 
event data for the different conversion mechanism 
as labelled in the figure. The error displayed in 
the experimental points is only statistical. 
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Figure 3. Angular distribution for Super- 
Kamiokande upgoing muon data together with 
our prediction in the absence of oscillation as well 
as the prediction for the best fit point to the full 
data sample for the different conversion mecha- 
nism as labelled in the figure. 
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Figure 4. Allowed regions of th neutrino decay 
parameters for the different Super-Kamiokande 
data samples. 



in the sub-GeV and multi-GeV samples present 
a similar deficit. As a consequence the allowed 
decay lifetimes which give a good description to 
the sub-GeV and multi-GeV data present very lit- 
tle overlap as can be observed by comparing the 
upper and central panels in Fig. (^. As a con- 
sequence the mechanism gives a worse fit to the 
global contained sample. This can also be ob- 
served in the angular distribution of contained 
events for the best fit to the contained events 
(dash-dotted line in Fig. (||)) from where one sees 
that the decay hypothesis cannot produce enough 
up-down asymmetry for the multi-GeV sample 
without conflicting with the sub-GeV data. 

This behaviour becomes particularly lethal 
when trying to describe the upward going muon 
data since for lifetimes favoured by the contained 
event data very little muon conversion is expected 
already for stopping muons in contradiction with 
observation. Based on this fact this mechanism 
was ruled out in its simpler form in Ref. |23 2J| 



Recently the possibility of neutrino decay in a 
more general four neutrino scenario has been re- 
visited in Ref. P6[ where it is discussed that a 
good description to the full atmospheric data 
sample is possible. 

I finally discuss our results on the alternative 
explanation of the atmospheric neutrino data in 
terms of FC neutrino-matter interactions iQ . In 
Fig. (^ I show the contours of the regions allowed 
by the Super-Kamiokande data. The different 
panels of the figure refer to the fits performed over 
the different sets of data. The shaded areas are 
the regions allowed at 90% C.L., while the dashed 
and dotted contours refer to 95 and 99 % C.L., 
respectively. In Table I list the corresponding 
best fit points as well as the values of Xmin attain- 
able for this mechanism. From the table we see 
that contained and stopping muon events can be 
described with a quality comparable to the oscil- 
lation channels. In Fig.(0) one can also observed 
that the FC-neutrino mechanism leads to a very 
good description of the zenith angular distribu- 
tions for contained events. This may come as a 
surprise since angular distributions for multi-GeV 
and sub-GeV events are rather different while the 
FC mechanism leads to an energy independent 
conversion probability. One must bear in mind, 



however, that the plotted angular distribution is 
that of the produced charged leptons in the neu- 
trino interaction. For neutrinos leading to sub- 
GeV events the average opening angle between 
the neutrino and the produced lepton is 60 de- 
grees what leads to the flat observed distribution 
almost independently of the specific conversion 
mechanism. The overall normalization, is how- 
ever, totally consistent with the FC hypothesis 
as the deficit in both samples is of about 60 %. 

The allowed regions can be qualitatively un- 
derstood in the approximation of constant matter 
density in Eq.(p^). From there one can see that 
in order to have a relatively large transition prob- 
ability, as implied by the contained events and 
also by the stopping muons events, the FC pa- 
rameters are required to be in the region ej, < Ci, 
and rj > Tr/i?®. This last condition leads to a 
lower bound on e^- The island in Fig. da.a) cor- 
responds to 77 ~ tt/Rq. On the other hand, the 
through-going muon data require a smaller tran- 
sition probability and therefore the region e'^ > e^ 
turns out to be the preferred one. 

The combination of the different data sets in 
a single x^-analysis is shown in Fig. (p[d). As 
seen in the figure as well as in Table ^ when 
the angular information of both stopping and 
through-going muons is included in the data anal- 
ysis, the quality of the full description worsens 
and leads to a X™,„ = 43.8/(33d.o./). Since FC- 
neutrino interactions lead to an energy indepen- 
dent conversion probability, the smaller deficit in 
the through-going muon sample as compared to 
contained and stopping-muon samples cannot be 
well accommodated. As seen in the lower panel 
in Fig. (0) the prediction for the best fit point 
for this mechanism would imply a larger deficit 
for higher energy muons mainly in the last three 
angular bins where it does not produce a suffi- 
cient amount of through-going muons at angles 
< 9 < 20 degrees below the horizon. 




Figure 5. Allowed regions of the FC-ncutrino 
interactions parameters for the different Super- 
Kamiokande data samples. 
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Table 1 

Minimum value of x^ and the best fit point for each channel and for different data sets. 



Experiment 


Oscillation 




jy^ -^ Vr 


Am'^ < 


I'll -^ Vs 

Am^ > 


FC i^-matter 


Super-Kam 






8.8 


12.9 


12.6 


9.4 


contained 


Am2 ( 10-3eV2 


) 


2.5 


3.2 


3.0 


6=0.95 


d.o.f=20-2 


sin^ 29 




0.98 


1. 


0.99 


e'=0.084 


Super-Kam 






1.3 


2.4 


2.3 


1 


Stopping-^ 


Am^ ( lO-^eV^ 


) 


3.0 


3.3 


3.7 


£=0.76 


d.o.f=5-2 


sin^ 29 




0.99 


1. 


0.93 


e'=0.19 


Super-Kam 


Am«n 




10.4 


13.4 


10.5 


10.4 


Through- Going-/i 


Am2 ( 10-3eV2 


) 


10.0 


4.9 


18. 


6=0.08 


d.o.f=10-2 


sin^ 29 




0.78 


1. 


0.55 


e'=0.26 


Super-Kam 


Amtn 




23.5 


32.9 


32.5 


43.8 


Combined 


Am2 ( 10^3gV2 


) 


3.1 


4.1 


4.5 


e =0.57 


d.o.f=35-2 


sin^ 29 




0.98 


1. 


0.96 


e'=0.45 



